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Here, we study noisy transitional flows in imperfect millimeter-scale channels. For probing 
the flows, we use microcantilever sensors embedded in the channel walls. We perform 
experiments in two nominally identical channels. The different set of imperfections in 
the two channels result in two random flows in which high-order moments of near¬ 
wall fluctuations differ by orders of magnitude. Surprisingly however, the lowest order 
statistics in both cases appear qualitatively similar and can be described by a proposed 
noisy Landau equation for a slow mode. The noise, regardless of its origin, regularizes 
the Landau singularity of the relaxation time and makes transitions driven by different 
noise sources appear similar. 


1. Introduction 

More than a century ago, Osborne Reynolds investigated transition to turbulence 
in a glass pipe in which he injected a filament of dye at the inlet (Reynolds 1883). 
Reynolds noticed that the characteristics of the dye filament and hence the entire flow 
field depended on the dimensionless flow rate or the Reynolds number, Re = UD/u (here, 
U is the mean flow velocity, D the pipe diameter, and v the fluid kinematic viscosity). 
When Re was below a critical value, Re Re cr , the dye propagated downstream as a 
thin filament without breaking up, indicating a laminar flow in the pipe. At Re ^ Re cr , 
this pattern changed dramatically: waves appeared in the vicinity of the inlet, sometimes 
leading to substantial flow randomization downstream. With increasing Re, the fraction 
of the tube length occupied by these waves increased, and at Re Re cr , the entire flow 
became turbulent. 

Reynolds, however, was unable to determine a unique value for Re cr . He noticed that 
the value of Re cr was sensitive to various imperfections, most notably the geometry of the 
inlet. If the inlet was sharp, inlet perturbations appeared in the form of shedded vortices 
that caused transition to turbulence at large Reynolds numbers. These perturbations, 
however, rapidly decayed downstream, if Reynolds number was not too large. By carefully 
eliminating these, Reynolds was able to delay the transition up to Re ~ 12, 800 (Reynolds 
1883). Others following Reynolds were able to sustain laminar flows in pipes even for 
Reynolds numbers as high as 100,000 (Pfenninger 1961). Relatively recently, the effects 
of initial (inlet) perturbations on Re cr were quantified by introducing well-controlled 
disturbances (jets) injected through holes in the vicinity of the inlet (Darbyshire & Mullin 
1995). Consistent with Reynolds’ observations, the flow became turbulent at smaller and 
smaller Reynolds numbers as the ratio of the disturbance velocity to the mean flow 
velocity was increased. (For a comprehensive review, see Yaglom (2012).) 
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The above-described “decay or amplification” of waves (or perturbations) forms the 
basis of Landau’s phenomenological theory of transition (Landau & Lifshitz 1987). It 
must be stressed that, based on his general approach which revolutionized the theory 
of critical phenomena, Landau was not interested in the notoriously difficult and non- 
universal problem of deriving the value of the critical (transitional) Reynolds number 
Re cr . Instead, he assumed the existence of a transitional (marginal) velocity distribution u 
at Re = Re cr , and attempted an investigation of the general and universal flow properties 
in the vicinity of Re cr in terms of a small perturbation velocity v. Then, under the small 
(infinitesimal) perturbation, the total velocity becomes u = u + v and the total pressure 
becomes p = po+Pi- I n general, the field u can be time-dependent but, following Landau, 
we assume the transitional pattern to be steady. The Navier-Stokes equations can then 
be written as 

u • Vu = + zW 2 u, V • u = 0; (1.1) 

P 


—^ + u • Vv + V • Vu +V • Vv =- — + v\7 2 v, V • v = 0. (1.2) 

at p 

Both u and v vanish on solid walls. Neglecting the 0{v 2 ) contribution to Eq. (1.2) 
results in linearized Navier-Stokes equations to be used for investigating instabilities 
in fluid flows. Now, the task is to find a solution to Eqs. (1.1-1.2) describing the time 
evolution of an initially (t = 0) infinitesimally small perturbation v. In this case, the 
0{v 2 ) contribution to Eq. (1.2) is neglected. For the perturbation, Landau assumed the 
form v(r, t) = A(t) f(r), where A(t) = constant x e~ li2t is the slowly-varying amplitude 
with a complex eigenfrequency, Q = uj + iy, and f(r) describes the spatial dependence. 
If the near-wall effects are not too important and in the frame of reference moving with 
mean flow, one can write the differential equation 


44J2 

dt 


= 2 7l |A| 2 -a|A| 4 , 


(1.3) 


which is to be solved subject to initial perturbation A(0) = Aq. The solution to this 
equation is 

\ml = e 27lt n 41 

I Ail 2 l+|A| 2 ^(e 2 ^-l)' K - } 

Indeed, setting 7 = c(Re—Re cr ), one can reproduce the observed “decay or amplification” 
of perturbations. When Re — R cr < 0, A(t) —» 0 in the limit t —» 00 (t I/ 7 ); otherwise, 
the amplitude grows, saturating at |A(oo)| oc >/Re — Re cr . The theory suggests that 
information about the initial conditions disappears on a time scale r = oc | . 

This behavior is similar to the “critical slowing down” in the proximity of a critical 
point in phase transitions. In pipe flows, it has important and interesting consequences. 
A perturbation occurring at position l and being advected by a mean flow of velocity U 
stays in the pipe for a time interval to ~ (L — l)/U , where L is the pipe length. Therefore, 
to observe the decay of a perturbation generated in the bulk or its growth toward a final 
turbulent state, one needs to r, requiring unreasonably long pipes around Re cr . The 
divergent relaxation time r oc |Re — R cr | - ^ with ( ~ 0.56 was also recently obtained in 
numerical simulations of transition in force-driven Navier-Stokes equations by McComb 
et al (2014). 

Landau’s theory of transition, though insightful, is better applicable to wakes in flows 
past bluff bodies (Sreenivasan et al. 1987) and in convection, i.e., in situations where 
wall effects and viscosity do not dominate. It is well-justified in pipe flows when the 
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characteristic size of a perturbation is O(D) and wall effects are unimportant. Efforts 
to describe transition in pipes using the Landau theory, most notably by Stuart and 
others (Stuart 1971), focused on obtaining the magnitudes of the coefficients 71 and a 
and testing their possible universality. This universality remains elusive, suggesting that 
wall effects must be important in transition. Numerical simulations and experimental 
data show that, at least in the range 2300 ^ Re ^ 10 5 , powerful bursts generated by 
unstable boundary layers are mainly responsible for turbulence production in the bulk. 

The majority of workers studying transition to turbulence in pipes have been interested 
in the flow response to perturbations in otherwise perfect pipes (Yaglom 2012). This inter¬ 
est can partially be explained by the mathematical well-posedness of the problem and by 
the emergence of numerical methods combined with powerful computers. Conversely, the 
“fuzzy” problem involving inlet disturbances, pipe imperfections, and pipe roughness has 
not attracted as much attention (Pausch & Eckhardt 2015). In this article, we investigate 
both experimentally and theoretically transition to turbulence in imperfect channels. In 
other words, we are not interested in reducing roughness or removing wall and inlet 
imperfections. Specifically, we strive to quantify growing perturbations near the wall by 
their spectra and statistical properties, including probability densities and low- and high- 
order moments. Remarkably, we observe that, while low-order moments are relatively 
independent of the experimental details, high-order moments are extremely sensitive 
to these details and may differ by orders of magnitudes. To describe our experimental 
observations, we propose a modified Landau theory in which all imperfection-induced 
flow disturbances are treated as added high-frequency noise. This theory agrees well 
with our experimental observations. An important consequence of this theory is that the 
noise regularizes the Landau singularity of the relaxation time. It must be emphasized 
that, in this regime (i.e., in imperfect channels), turbulence is constantly driven by noise, 
which makes the problem very different from the initial value problem considered by Hof 
et al. (2006). 


2. Experiments 

2 . 1 . Experimental Set-up 

We have performed our experiments in two nominally identical rectangular channels 
with linear dimensions L x x L y x L z = 35 x 8 x 0.9 mm 3 and with hydraulic diameters 
D & 1.6 mm. A pressurized air source is used to establish air flow in each channel. 
The flow rate is monitored using a rotameter, and the pressure drop Ap between the 
inlet and the outlet is measured using a differential gauge. The near-wall flow in each 
channel is probed by a rectangular microcantilever with linear dimensions l x x l y x l z = 
27.5 x 225 x 3 /am 3 [Fig. l(a)-(b)]. In the first channel [henceforth channel (i)], the chip 
carrying the microcantilever is embedded in the bottom wall of the micro-channel by 
a surrounding aluminium structure that is cut from a ~ 300 — /im-thick smooth sheet 
(matching the height of the chip) and glued to the bottom wall [Fig. 1(b)]. The root 
mean square (r.m.s.) roughness on the wall is ^ 300 nm. In the second [henceforth 
channel (ii)], a groove, which matches the size of the cantilever chip closely, is machined 
on the bottom channel wall [Fig. 1(b)]. The r.m.s. roughness on this wall is significantly 
larger, ~ 3 /am. In both channels, the test section is approximately 17 mm (« 11 D) from 
the inlet. Because turbulence is driven and sustained by noise due to imperfections in 
our channels, the distance of the test section to the inlet is sufficient for the observation 
of the transition (Zagarola & Smits 1998). Both channels have smooth and transparent 
top walls, and the motion of the tips of the microcantilevers is read out using the optical 



Figure 1. (a) Rectangular channel with linear dimensions L x x L y x L z = 35 x 8 x 0.9 mm 3 . The 
test section is highlighted. The microcantilever chip is integrated to the middle of the bottom 
wall, (b) Optical microscope (left) and scanning white light interferometry images of the bottom 
walls of the two channels. The optical images are 6.2 x 4.6 mm 2 , and the interferometry images 
are 0.7 x 0.5 mm 2 . The color scale bar is the same for both images, (c) Double logarithmic plot 
of pressure drop Ap as a function of Reynolds number in the channels. Lines are fits to laminar 
(dashed) and Blasius (dotted) flow models with the transition around Re ^ 2 x 10 3 . Error bars 
are shown only when larger than symbols. 


beam deflection technique (Meyer & Amer 1988). We confirm that our optical transducer 
remains linear in the explored parameter space. 

We emphasize that the two channels, while nearly identical in all macroscopic aspects 
(e.g., linear dimensions, inlet pipes), possess different sets of microscopic imperfections 
on their walls and test sections, as seen in Fig. 1(b). Here, we are not concerned with 
the particulars of these imperfections, which may include but are not limited to surface 
roughness, wall asperities, bumps, edges and so on. As we show below, even though these 
imperfections may lead to different random flows, the lowest order statistics remain 
qualitatively similar. 


2.2. Pressure Drop 

The pressure drop Ap between the inlets and outlets of both channels are shown in 
Fig. 1(c) as a function of the Reynolds number. In the low Reynolds number range 
0 ^ Re < 10 3 , Ap can be fit to that in a Poiseuille flow using the nominal channel 
parameters with no free parameters. This is the dashed asympote in Fig. 1(c). In the high 
Re regime, the data appears to converge to another asymptote (dotted line). This is the 
pressure drop calculated from the Colebrook equation in a rectangular duct (Jones 1976), 
using only experimentally available parameters. The pressure drops in both channels 
match closely. 


2.3. Spectral Properties 

We first turn to the spectral properties of near-wall fluctuations. Flow forces act on 
the microcantilever and give rise to mechanical fluctuations. The microcantilever has a 
sharply-peaked resonance at 100 kHz with a linewidth of 500 Hz. Its linear response 
function is frequency-independent in the frequency range / < 80 kHz and can be 
approximated as |G(/)| ~ 1/ft, with k being the cantilever spring constant and k ~ 3 
N/m. Thus, the linear relation, S z (f ) ~ Sf(/)/ ft 2 , between the spectral densities of the 
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Figure 2. Power spectral density S z (f) of cantilever displacements for different Reynolds 
numbers for channels (i) and (ii). The dotted line in (a) shows the noise floor of the measurement; 
the noise floor and the spectra are indistinguishable at Re ~ 820 for the frequency range / ^ 300 
Hz. 


force Sf (/) and the cantilever displacement S z (f) remains valid below / < 80 kHz. At 
very low flow rates (Re < 800), S z (f) becomes obscured by measurement noise because 
the flow cannot generate detectable cantilever motion. 

The spectral densities of the cantilever fluctuations S z (f) obtained for the two channels 
are shown at different Reynolds numbers in Fig. 2. For Re < 1000 [Fig. 2(a) and (b)], 
the spectra for both cases are barely above the noise floor (dotted black line) and appear 
similar for the most part. For Re ^ 1600 [Fig. 2(c)-(f)], small differences between the 
two cases can be noticed. First, the fluctuations in channel (ii) with the rougher wall 
are larger than those in channel (i) by an order of magnitude; second, the spectrum 
in channel (ii) extends to slightly higher frequencies. As Reynolds number is increased, 
both data traces increase monotonically and smoothly. Interestingly, S z (f ) oc 1 // for 
both channels in the low frequency range. 
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Normalized Displacement z/z rms 

Figure 3. PDFs at different Reynolds numbers for channels (i) and (ii). The x-axes are in 
units of Zrms , which are different for each data set (see table 1). Dotted lines are Gaussians. 


2.4. Statistical Properties 

Next, we turn to the statistical properties of the near-wall fluctuations. We collect 
long-time traces of the cantilever amplitudes filtered in the frequency range 100 Hz < 
/ < 30 kHz to remove the high-frequency resonant oscillations. We then sample the 
time data every 3 //s, obtaining ~ 10 6 data points. Because we are interested in single¬ 
point probability density functions (PDFs), we do not worry about over-sampling in 
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comparison to other relaxation times in the measurement, e.g., the ring-down time of the 
cantilever. 

The PDFs for different Reynolds numbers are plotted in Fig. 3. For Re < 1000, our 
measurements are dominated by technical noise, and the PDFs in both channels (i) and 
(ii) are perfectly Gaussian [Fig. 3(a) and (b)]. In the range Re ^ 1600, both PDFs 
go through dramatic changes as seen in Fig. 3(c)-(f). First, one notices a substantial 
broadening of the tails of the PDF in channel (i), indicating the presence of strong wall 
velocity/pressure bursts, which may reach the bulk (Yakhot et al. 2010). The PDFs 
observed in channel (ii) for Re > 1600 suggest that the flow here is somewhat more 
homogenous [Fig. 3(c)-(f)] compared to that in channel (i), but wall bursts make the 
PDFs more asymmetric and intermittent [Fig. 3(e)]. This trend continues with increasing 
Reynolds number [Fig. 3(f)]. The dotted lines in the insets show that the PDFs in channel 

(i) cannot be fit to Gaussians, even at very small displacements; an exponential decay 
seems perhaps more appropriate. The differences between the two flows can also be clearly 
seen in the moments of the PDFs listed in Table 1. The moments are also plotted in Fig. 4, 
with z rms = (z 2 ) (brackets indicate averaging) in Fig. 4 (a) and normalized high-order 
moments in Fig. 4 (b). The high-order moments obtained in the two channels differ by 
orders of magnitude. These observations can be summarized as follows. The random 
flow in channel (i) is more intermittent compared to that in channel (ii); however, the 
intensity of fluctuations in channel (ii) is larger. Based upon these observations, it may 
not be incorrect to conclude that these are two very different random flows. While the 
exact source of the differences in these two flows is impossible to identify, it must be due 
to the different set of imperfections in the two channels. 

We return to the r.m.s displacement of the cantilever, z rrns = (z 2 ) , as a function of 

Reynolds number, shown in Fig. 4 (a). The dashed lines are fits to our proposed equation, 
discussed below. We note that z rrns values obtained by integrating the spectra and from 
the time domain measurements agree closely. There are three regions in the plot marked 
with different shadings. The first region, Re < 800, is dominated by technical noise and 
does not provide any insight. In the second region (blue) where 1000 < Re < 1600, the 
magnitude of the cantilever fluctuations are of the same order for both channels (i) and 

(ii) , z rrns ~ 1 nm. In the third region (pink), the observed r.m.s. amplitude of near-wall 
fluctuations in the rough channel (ii) are larger than those in the smooth channel (i). 
The slope changes for both data traces around Re « 2000, suggesting the onset of more 
sustained perturbations. The data traces appear to increase parallel to each other for 
Re > 2000. 


3. Theory 

The similarities in the two data traces in Fig. 4(a) suggest that there may be common 
underlying physics to both cases. In order to describe both flows by a single equation, we 
return to the Landau theory discussed above. We incorporate the non-idealities into the 
Landau theory by considering a general additive noise term </>(£). We assume that a high- 
frequency random Gaussian force </>(£) defined by the correlation function </>(£)(/>(£') = 
2 <j) 2 5{t — t') stirs the fluid. Then, Eq. (1.2) is modified as 


— + u • Vv + v • Vu + v • Vv = 

at 


Vgi + 


zW 2 v + 


P 


(3.1) 
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15.5 
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0.07 

0.07 

3.00 

3.04 

15.0 

15.4 

105 
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820 

0.07 

0.08 

2.99 

3.01 

15.0 

15.0 

106 

103 

1030 

0.08 

0.07 

3.00 

3.06 

15.1 

15.6 

107 

111 

1240 

- 

0.14 

- 

3.75 

- 

29.2 

- 

366 

1650 

0.25 

0.56 

5.61 

3.92 

81.5 

31.8 

1890 

421 

2100 

0.97 

2.50 

8.75 

4.15 

261 

47.3 

15400 

1500 

2330 

1.13 

3.33 

13.7 

3.94 

654 

38.1 

60400 

768 

2550 

1.46 

4.05 

20.7 

4.34 

2840 

49.6 

1.02 x 10 6 

1120 

2750 

1.79 

5.02 

31.2 

4.65 

5960 

55.4 

2.36 x 10 6 

1170 

3000 

2.82 

5.73 

48.4 

5.33 

20100 

97.6 

1.39 x 10 7 

5010 


Table 1. Moments. 
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Figure 4. (a) Cantilever r.m.s. displacements, z rms , versus Re for channels (i) and (ii). Dashed 
lines are fits to the results of the noisy Landau Equation, Eq. (3.5), with parameters given in 
Table 2. (b) Normalized high-order moments. Error bars are estimated from convergence tests 
performed by computing the moments for a subset of the data (e.g., half). When not shown, 
error bars are smaller than symbol size. 


Repeating Landau’s arguments leading to Eq. (1.3), we write 


d A . a . o ,, x 

-= 1iA --a +m 


(3.2) 


Averaging over high-frequency fluctuations gives a modified Landau equation for the slow 
mode 

d\A\ 2 


dt 


= 2 7l |A| 2 -a|A| 4 + 0 2 . 


(3-3) 


In a force-driven flow with initial condition Aq = 0, the solution to Eq. (3.3) is 


\A(t)f = 21 + p - ~ P +e ~ 
a /?- + /?+ e ” 


(3.4) 
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Re cr 71 /a (m 2 /s 2 ) (j) 2 /a (m 4 /s 4 ) 


Channel (i) 1850 0.14 x (Re - Re cr ) 2.37 x 10 3 
Channel (ii) 1750 0.29 x (Re - Re cr ) 5.41 x 10 3 

Table 2. Fit parameters used [see Eq. (3.5)]. 


where fl± = d= ^ + \ ^ + — > 0 with r = 

i -L- rv V oi 


2 V ll+<t> 2 u 


This gives in the long-time limit 


\A(oo)\ 2 — — — + \l + —. 


/ 71 2 (j) 2 


a 


a* a 


(3.5) 


Remembering that in Landau theory 71 = c(Re — Re cr ), we notice an important conse¬ 
quence of noise: the relaxation time r remains finite in the limit Re —>> Re cr , in contrast to 
the “critical slowing down” discussed above. Thus, the external noise source regularizes 
the dynamics around the transition point. In addition, Eq. (3.2) indicates that, when A 
is small in the low Re limit so that the G(A 3 ) contributions can be neglected, A obeys 
Gaussian statistics. 

The obtained results can be used to quantitatively explain the experimental data of 
Fig. 4(a). The time-dependent force acting on the cantilever is 

F = Fz = - [ pi • n dS& T-Vz, (3.6) 

Js dz 


where z is the unit vector, and V is the volume of the cantilever. We simplify the 
problem by approximating our channel as a long and wide rectangular (planar) duct 
and, neglecting noise, write the perturbation equation in the z direction from Eq. (3.1) 
as 


1 dpi 
p dz 


dv z 


dt 


4z 2 \ dv 


dv 2 N 


St'VV + OIx 


dz 


(3.7) 


In close proximity to the wall, z ~ H/ 2 , and the second term on the right hand side of 
Eq. (3.7) is small. Numerical simulations (Lee et al. 2004), where statistics of acceleration 
in close proximity to a wall was studied, suggest that bursts are dominated by the z- 
component of the velocity ( v 2 = 0(v z 2 )) and the viscous term is unimportant. Therefore, 


l^Pi _dv L 

p dz dt dz 


(3.8) 


Here, d t v z = 0(v z U / H). To find an order of magnitude estimate for d z {y 2 \ we 
extrapolate the results of Yakhot et al. (2010), which shows that, for Re < 10 5 , the 
velocity of wall bursts are v z = 0(77/10). This means that the magnitude of the two 
terms on the right hand side of Eq. (3.8) are comparable, and 0(v z 2 ) must be accounted 
for. Thus, we deduce d z p\ = 0{v z 2 /H ) and write 


Zrms ~ pVv z 2 /kH OC |H(oc)| 2 . 


(3.9) 


This exercise provides the fits (dashed lines) shown in Fig. 4 (a) with the parameters in 
Table 2. The ultimate justification for the above arguments comes from the agreement 
between experiment and theory in Fig. 4 (a). The following simple order of magnitude 
estimate further bolsters our confidence. Assuming v 2 jU 2 ~ 0.008 around Re cr (Yakhot 
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et al. 2010) and using the experimentally available parameters (V ~ 2 x 10“ 14 m 3 , 
H ~ 10 -3 m, ft ~ 3 N/m), we find z rms ~ 10“ 10 m, close to the experimental values. 


4. Summary and Outlook 

In transitional flows in imperfect channels, high-order moments of near-wall fluctua¬ 
tions appear to be extremely sensitive to the specifics of the imperfections in the channels. 
In the two nominally identical channels with different sets of imperfections presented 
here, the high-order moments of the near-wall fluctuations differ by orders of magnitude. 
On the other hand, low-order moments remain relatively independent of the details of 
the imperfections. In our experiments, z rms values seem to depend directly on the wall 
roughness. Additionally, z rrns in these flows can be accurately described by the noisy 
Landau equation presented. The applicability of the noisy Landau equation to these 
different flows is probably due to the fact that the noise term is taken to be at high 
frequencies as compared to the slow mode in question. In other words, the stirring of 
the fluid occurs at high frequencies. This assumption is eminently reasonable because 
microscopic surface asperities giving rise to roughness are at high spatial frequency 
compared to any length scales of the flow. Our preliminary results suggest that inlet 
disturbances may also be accounted for by the noisy Landau equation. Along the same 
lines, we believe that even thermal fluctuations in the fluid would likely result in a similar 
regularization of the Landau equation. Finally, the phenomena observed here may have 
important consequences for heat and mass transfer in wall-bounded flows; these will be 
discussed in detail in future work. 

We acknowledge support from the US NSF through grant Nos. CMMI-0970071 and 
DGE-1247312. We are grateful to T. Kouh for intial measurements, and K. R. Sreenivasan 
and J. Schumacher for valuable comments. 
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